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Ab?tra<?£ 

The paper presents the 
verification of a new method for 
computing sensitivity derivatives of a 
coupled system. The method deals with a 
system whose analysis can be partitioned 
into subsets that correspond to 
disciplines and/or physical subsystems 
that exchange input-output data with each 
other. The method uses the partial 
sensitivity derivatives of the output 
with respect to input obtained for each 
subset separately to assemble a set of 
linear, simultaneous, algebraic equations 
that are solved for the derivatives of 
the coupled system response. This 
sensitivity analysis is verified using an 
example of a cantilever beam augmented 
with an active control system to limit 
the beam's dynamic displacements under an 
excitation force. The verification shows 
good agreement of the method with 
reference data obtained by a finite 
difference technique involving entire 
system analysis. The paper also 
demonstrates the usefulness of a new 
system sensitivity method in optimization 
applications by employing a piecewise- 
1 inear approach to the same numerical 
example. The new method's principal 
merits are its intrinsically superior 
accuracy in comparison with the finite 
difference technique, and its 
compatibility with the traditional 
division of work in complex engineering 
tasks among specialty groups. 

Ug* of Symbol? 

Ai vector of functions 

b beam rectangular cross-section 

width 

c damping constant - design variable 

vector for the controls subsystem; 
and, controls subsystem subscript 
f functional relationship 

F(t) dynamic forcing function 

G gain matrix 

h beam rectangular cross-section 

depth 

I identity matrix 

J(x,y) Jacobian matrix of dx/dy 
M mass (inertia) matrix 

m mass 

N number of the function vectors A^ 

ni length of vectors A^ and Yi 
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p Ricatti matrix 

r vector of time-dependent 

displacements 

s structures subsystem subscript 

Ti static transverse load on beam 

u control input variable 

w vector of static displacements 

x state vector 

X vector of design variables 

Y vector of unknowns in the system 

analysis 

rj generalized coordinate vector 

* matrix of natural vibration modes 

p user controlled factor in the 

Kreisselmeier-Steinhauser function 
a) natural vibration frequency 

critical damping for i-th 
vibration mode 

Other symbols used locally are defined 
whore first introduced. 

Introduction 

The merger of structures and active 
controls in control -augmented structures 
appears to be one of the most promising 
means that have recently become available 
to engineers for making a quantum jump in 
structural efficiency, especially for 
dynamic applications in flight and ground 
vehicles, and in space structures. While 
extending the design freedom, inclusion 
of the active controls confronts 
structures designers with a coupled 
system whose behavior is a resultant of 
the structures-controls interaction 
governed by the design variables 
available in both disciplines, 
complexity of that behavior limits the 
effectiveness of intuition and past 
experience as design guides and suggests 
the use of a formal sensitivity analysis 
to support human judgment and to provide 
a basis for numerical optimization. 

Recognizing the above, ref.l 
presents a method for sensitivity 
analysis of the control -augmented 
structure, and ref. 2 reports on the 
numerical optimization results using the 
sensitivity data. The sensitivity 
analysis developed in ref.l is based on 
the classical, quasi-analytical approach 
of representing the entire system by a 
single set of governing equations that 
incorporates terms pertaining to both 
subsystems - structures and controls. 
These equations are differentiated with 
respect to the design variables to yield 
the sensitivity equations containing 
the derivatives of behavior as unknowns. 
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The purpose of the study reported 
herein is to test a new, alternative 
method of evaluating the structure- 
control system behavior sensitivity 
derivatives using the decomposition 
approach introduced in ref. 3, and to show 
an example of the use of the method in 
optimization. In the method of ref. 3, 
the sensitivity equations are assembled 
from building blocks in which each block 
represents the partial sensitivity 
derivatives of each subsystem's output 
with respect to its input and its design 
variables. The principal merit of this 
alternative method is that the partial 
sensitivity data are obtained from 
separate subtasks, self-contained within 
each of the subsystems. This supports 
the division of large engineering 
projects among specialty groups, is 
compatible with the technology of 
distributed computing, and enables one to 
use experimental sensitivity data. 

Structures-Controls System 

It is expedient to use a test case 
to introduce the method. The test case 
structure is a simple, cantilever beam 
shown in Fig.l, subjected to static loads 
and to a dynamic excitation force at the 
tip. The model is controlled by two sets 
of actuators - one for the tip lateral 
displacement and another for the tip 
rotation. The active controllers limit 
the beam's dynamic displacements within 
prescribed bounds. The structure must 
also be sized to limit the static 
stresses below allowable levels. 

The optimization problem to be 
solved for this system calls for the 
minimization of two objectives: the sum 
of the structural weight of the beam and 
the weight of the control system, and the 
control system effort. The control 
system weight component is assumed to be 
a simple analytical function of the 
control effort. Minimization of the 
control system effort is carried out 
within the control subsystem by recourse 
to the classical linear quadratic 
controller synthesis. 

Constraints are imposed on the 
static displacements and stresses, 
dynamic displacements, and natural 
frequencies. The design variables are 
the ratio of critical damping to 
frequency in the control subsystem and 
the beam cross-sectional dimensions 
in the structure subsystem. Additional 
details of the example are available in 
the Appendix. 

In order to perform the 
optimization, it is necessary to 
calculate derivatives of the entire 
system response with respect to the 
design variables? the new method for 
sensitivity analysis is implemented for 
that calculation. 


Structures Subsystem 

The beam constitutes a structure 
subsystem in the structure-control 
system. It is discretized in the 
span-wise direction into five segments of 
equal length. In the finite element 
displacement method used in the study, 
the finite element model has six nodal 
points, twelve elastic degrees of freedom 
(two per nodal point) , and five two- 
dimensional beam finite elements, A 
lumped mass representation is used to 
model the beam inertial characteristics. 

The structural analysis for the beam 
dynamic response is comprised of 
equations for eigenfrequencies and 
eigenmodes 

(K - u> 2 M)*i = 0 (1) 

and of load-deflection equations for the 
static loads 

K w = T (2) 

where w is a vector of static 
displacements and T is the load vector. 
Two different sets of design variables 
are used in the structure subsystem. The 
first set includes the reciprocal of 
the cross-sectional bending moment of 
inertia and the cross-sectional area of 
each segment of the beam. The other 
combination is the width and depth of the 
rectangular beam cross-section for each 
segment. 

Controls Subsystem 

It is assumed that the actuators are 
part of an optimal control system with a 
linear quadratic regulator, so that the 
control system is described by the 
following equations. The equations of 
motion for the dynamical system are 

[ M ] { r } + [C]{f> + [K] { r}_ = 

[ b ] { u } + [b]{F} (3) 

A modal transformation, that can be 
written as 

{r> = [*){*) (4) 

enables one to define a state vector as 



in terms of which the state equation 
becomes 

(x) = [A] { x ) + [B]{u> + [5] <F) (6) 

where 


o 

H 


0 

L. 

[B] = 


L-U) 2 1 -2C w 

4> T b 
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Generic Solution 


[B] - 


0 

# T b 


(7) 


The optimal control problem is one 
in which the following quadratic 
performance index is minimized 

J - ({x) T [Q]{x)+(u}' r [R](u))dt (8) 

Jo 

subject to 

(X) - [A] (X) + [B]{u> 

where [Q] and [R] are arbitrary weighting 
matrices. The solution of the above 
linear quadratic regulator problem 
results in the nonlinear algebraic 
Riccati equation. 

[A] T [P] - [PHBHRr^B^tP] + 

[P][A] + [Q] ■ 0 (9) 


The solution of this equation yields the 
symmetric positive definite matrix 
referred to as the Riccati matrix , which 
is then used to obtain the gain matrix 
from the relationship 

[G] - IR]“ 1 [B] T [P] (10) 


Taking the approach from ref. 3, the 
system's governing equations are 
partitioned, with each partition 
corresponding to a distinct engineering 
discipline, and/or a physically distinct 
subsystem. The following notation 
defines the i-th partition 

= 0 (15) 

where Ai is a vector of functions equated 
to zero, so that they form a set of n^ 
equations. The quantities in the inner 
parentheses are the inputs that must be 
given in order to solve eq.15 for Y^. 
These inputs entail the independent 
variables X which include design 
variables and externally prescribed 
constants, and dependent variables Y j , 
j / i. The latter are obtained from 
solutions of the other partitions Aj , 
j / i, and represent couplings of Ai to 
Aj's. The couplings require that the 
entire set of governing equations 
composed of the partitions such as 
eq.15 must be solved as a set of 
simultaneous equations - an iterative 
solution method is commonly used for 
large, complex, and nonlinear systems. 


The control input vector, computed from 
the system response using the gain 
matrix, can now be written as 

(u) - - [G] ( x) (11) 

and substituted into the state equation 
so that the controlled system response 
may be obtained from 

{x> - [A] (x) + [B] { F) (12) 

where 

[A] = [A] - [B)[G] (13) 

Once the state vector has been found from 
integration of eq.12, the dynamic 
displacement vector may be obtained from 
eq.4. 

The design variable in the control 
system is the damping constant, c, 
defined as 

2 Ti 

C - (14) 

w i 


When the solution Y, composed of all 
the partitions Yi, is obtained for a 
given vector X, the sensitivity 
derivatives of Yi with respect to X are 
computed as a solution of a set of 
simultaneous, linear, algebraic equations 


I 


• 

• 


: 1 

• 


dYi 


-H 

4-1 

I • • Jij 


— 

► = « 

— 

• • • 


ax k 

• 


ax k 

• 

• • • 


dYj 


££j 

Jji . .1 





_ 


. dX k . 


. *x k . 


(16) 

In these equations, the matrix of 
coefficients has unities on the 
diagonal and its off-diagonal submatrices 
are the Jacobian matrices of the partial 
derivatives of output with respect to 
input for the coupled Ai's. Specifically, 
expressing Yj[ as an implicit function of 
the inputs to A^ 


for the i-th eigenmode. In the above 
analysis, this variable affects the state 
solution via the matrix [A] . 


Yi — f i (X, . . . , Yj , • . . ) , j — i (17) 
the Jacobian J-^j is 


SYstem.. Sensitj.vjt:y Analysi s 

Solution to the structure-control 
system is based on the algorithm 
introduced in ref. 3. The essence of the 
referenced approach is stated first, 
followed by its adaptation to the problem 
at hand. 


i j 



n^ x n j 


( 18 ) 


It consists of the partial sensitivity 
information, obtainable locally for each 
Yf by treating it as a function of given 
arguments according to eg. 17. 
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The right hand side vector in eq.16 
contains the partial derivatives of 
Yi with respect to a particular X^, one 
at a time, also computable locally from 
eq. 17. 

Application to Numerical Example 

Applying the approach described in 
ref. 3 to the example case, we define the 
structure-control system by a schematic 
representation shown in Fig. 2, where S 
represents the structures subsystem and C 
represents the controls subsystem. The 
vector Y includes interaction of one 
subsystem with the other. The physical 
influence of the controls subsystem on 
the structures subsystem is through the 
actuator forces. However, the system 
equations of motion (eq.3) and evaluation 
of the dynamic displacement constraints 
have been included in the analysis of the 
Control Subsystem. Consequently, in the 
partitioned mathematical model of the 
system, the controls-to-structures data 
channel needs to transmit only the 
control system mass data that has to be 
accounted for in the vibration analysis 
in order to compute natural vibration 
frequencies and modes. 

The other data channel represents 
the structure's influence on the control 
system through the vibration frequency 
and mode shapes. In addition to the data 
coupling the two subsystems, there is 
data output to the outside representing 
response of the entire system. This data 
is the output from each subsystem as 
follows: 

Structures: weight, natural 

frequencies and the values 
of the static stress and 
displacement constraints. 

Controls: weights and the values of 

the dynamic displacement 
constraints. 

Table 1 summarizes the coupling data and 
the design variables for both subsystems, 
and additional numerical details are 
given in the Appendix. 

The generic governing equations may 
be written, in this particular case, as 
partitioned into the structural analysis 
S and control analysis C, and coupled by 
the presence of the output of one in the 
input of the other: 

S((X s ,Y c ) ,Y S ) = 0 (19) 

C((X c ,Y s ) ,Y C ) = 0 (20) 

These equations correspond to the groups 
of equations 1-2 and 3 through 13, 
respectively. 

When taken separately, equations 19 
and 20 yield the interactions, Y, as 
implicit functions of the design 
variables X s and X c as 


Y s * fs(Xs,¥ c ) ' Y c - f c (X c ,Y s ) ( 21 ) 

When solved simultaneously, equations 19 
and 20 produce the Y unknowns as 
functions of the variables X s and X c . 

Accounting for eq. 21, the generic 
sensitivities of eq.16 take on the 
following form: 

( 22 ) 


1 1 


dY s ] 


af* 

I |-J(f s ,Y c ) 


— 


— 

1 


dX* 


ax* 

1 


3Y C 



-J(f c ,Y s )\ I 


— 


0 

\ 


dX* , 




where the character 1 * 1 may be 
substituted by either X s or X c . 

The sensitivity solution, dY/dX, is 
obtained for one design variable at a 
time. 

The partial derivative data that 
make up eq.22 is obtainable 
separately within each of the 
participating "black boxes" of 
structures and control by any technique 
available, e.g., analytically, by finite 
differences, or could even be produced 
experimentally. A simple, one-step- 
forward finite difference technique is 
used in this study in both the structures 
and controls subsystems. 

Verification of the System Sensitivity 
Analysis and Results 

To verify the system sensitivity 
method represented by eq.22, the 
following stepwise procedure was 
implemented: 

1. The structure-control system is 
analyzed by solving eq. 19, 20 
for Y's for assumed values of X. 

2. Partial sensitivity derivatives with 
respect to X 1 s and Y 1 s are computed by 
finite difference for structures and 
controls separately, using, as 
reference values, the X's and Y's 
satisfying equations 19 and 20. 

3. The partial sensitivity derivatives 
obtained above are input to the system 
sensitivity eq.22, which are then 
solved for the sensitivity derivatives 
with respect to the X's. 

4 . The system sensitivity derivatives of 
Y with respect to X were calculated 
again by a finite difference 
technique, this time involving the 
coupled system analysis, i.e., eq.19, 
20 and using the same Y and X 
reference values as in step 2 . 

5. The system sensitivity derivatives 
obtained in steps 3 and 4 are 
compared. 
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This comparison is illustrated by 
typical data samples in Tables 2, 3, and 
4 that display sensitivity derivatives 
obtained for the following data: 

Beam length 500 in 

Beam Moment of Inertia 144 in4 

Beam Cross-sectional Area 48 in2 
Damping constant, c 0.01 

material properties for Al-alloy 

The constant parameters (gain 
levels) are set to achieve very low 
(almost negligible), moderate, and high 
levels of the control presence in the 
structural dynamics by changing the 
weighting matrices in the controls 
performance index. The results are shown 
in Tables 2, 3, and 4, respectively. 

All three tables pertain to the 
controls design variable and the columns 
in each table are, from the left: results 
from Step 4, results from Step 3, and 
relative error of the results from Step 3 
with respect to the results from Step 4. 

The sample shows that the agreement 
of the results produced by the method 
proposed in ref. 3, represented herein by 
eq.22, agree very well with the finite 
difference verification. The only 
discrepancies correspond to the 
relatively small derivative values, and 
are attributed not to error in the method 
being tested, but to a loss of accuracy 
of the finite difference technique caused 
by the limited word length in the 
computer. The sample shows also that the 
errors relative to the finite difference 
results are greater for the case of 
moderate level of control (Table 3) . 

This is so because the effect of the 
small finite difference step gets 
submerged in the numerical error of 
the simultaneous solution of eq.19, 20 - 
a confirmation of the effect predicted in 
ref. 3. The comparisons for other design 
variables of the problem were found to be 
qualitatively the same. 

gygtem Sensit ivity . Derivative? 

Used in , Ppt jin i g 3 1 i QT 1 

System sensitivity derivatives 
quantify "what if" questions that are an 
intrinsic part of the design process. 

They can also be used to guide a formal 
optimization algorithm. To demonstrate 
their usefulness in the latter, the 
structure-controls system was optimized 
by a piecewise-linear optimization 
method, using the system sensitivity 
derivatives calculated via eq.22 at each 
linear stage. 

The mathematical statement of the 
optimization problem is given in the 
Appendix, and the optimization procedure 
follows a sequence in which the system is 
initialized and a solution of this 
initial system (eq.19, 20) is obtained. 
The partial sensitivity derivatives for 
the structure and control subsystems 
are computed separately (see Step 2 in 


the procedure listed in the previous 
section) , followed by a solution of eq. 

22 for derivatives of the system response 
with respect to the design variables. 

The system is then optimized by 
approximating the objective function and 
constraints by linear extrapolation with 
respect to design variables, using the 
derivatives obtained above (see the 
Appendix for details of the formulation) , 
and restricting the change of design 
variables by prescribing lower and upper 
bounds. The results from this design are 
used as the new set of variables for 
which the sensitivities must be obtained. 
This procedure is repeated until 
constraints are satisfied and there is no 
appreciable change in the objective 
function for three cycles of design. 

The optimization procedure performed 
as expected. Fig. 3 and 4 illustrate 
typical histograms of the objective 
function of weight vs. iteration number 
for a feasible and infeasible initial 
design point, respectively. The usable- 
feasible directions algorithm (ref. 4) is 
used for searching the design space in 
Step 5. The overshoot peak of the 
histogram for optimization started from 
an infeasible point (Fig. 4) is typical 
for the usable-feasible directions 
algorithm. Although the problem is highly 
nonlinear, the procedure converged to 
essentially the same final design (see 
Appendix) confirming the system 
sensitivity analysis formulated in ref. 3 
as useful in optimization. 

Conclusions 

A method proposed in ref. 3 for 
computing sensitivity derivatives of an 
internally coupled system to design 
variables was tested by application to a 
simple, control-augmented structure. The 
method yields the system sensitivity from 
the solution of linear, algebraic 
equations built from the partial 
sensitivity information obtained 
independently for each of the subsystems. 

Comparison of the results with the 
reference results produced by a finite 
difference technique applied to the 
solution of the entire system confirmed 
good performance of the method of ref. 3. 
Also, it demonstrated that the finite 
difference technique when applied to an 
internally coupled system may lose 
accuracy due to the system solution noise 
overwhelming the effect of the small 
finite difference step. 

The usefulness of the sensitivity 
derivatives obtained in a manner 
presented in the paper was demonstrated 
by a successful application in the formal 
optimization of a control -augmented 
structure. The principal benefits of the 
system sensitivity analysis implemented 
in this paper are an improved accuracy in 
comparison to the finite difference 
method involving the entire system 
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analysis, and further, the decomposition 
of the system sensitivity analysis into 
partial sensitivity analyses performed 
separately in each discipline and 
subsystem. This should support a broad 
workfront of people and computers in 
design organization and allow use of 
specialized sensitivity analysis methods 
unique to each discipline. 


The following extrapolations are made in 
approximate analysis coupled to an 
optimization program: 


w t i+1 


. d(w s + w c )i 

Wfc 1 + AX S + 

dX s 
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Appendix 

Definition of the Example 

The optimization problem is stated 
minimize 

W t (X) = W S (X) + W C (X) (Al) 

subject to 

(A2 ) 

g(X) = f (static tip displacements) 
g(X) = f (natural frequencies) 
g(X) = f (static stresses) 
g(X) = f (dynamic displacements) 

X L < X < X 0 (A3) 


Move limits ranging from 2.5 % to 10 % 
are included in eq. A3. The constraints 
represented in eq. A2 are as follows: 

dl 

91 = - 1 < 0 (A6 ) 

dl a 

dr 

92 = - 1 < 0 

dr a 

“1 

g 3 = - + 1 < 0 

“la 

u>2 

g4 + i < o 

“2a 


95 = - 1 < 0 

a a 

ddl 

96 = - 1 < 0 

ddl a 

ddr 

g7 = — 1 < 0 

ddr a 

where the limits are set as: 


dl 

static lateral 
displacement 

50.0 

( A7 ) 
in 

dr 

static rotational 
displacement 

0.2 

rad 

“1 

first natural 
frequency 

1.0 

Hz 

a>2 

second natural 
frequency 

1.25 

Hz 

a 

static stress 

30,000 

psi 

ddl 

dynamic lateral 
displacement 

50.0 

in 

ddr 

dynamic rotational 
displacement 

0.15 

rad 


The stress constraint gs is formulated as 
a single cumulative constraint 
representing static stress in all the 
elements of the beam due to static loads 
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PPtiJtigation Results 


Ti, T 2 , T 3 . The cumulative constraint 
is taken in form of the Kreisselmeier- 
Steinhauser function (KS function) 
introduced in ref. 5. The KS is a 
continuous, differentiable envelope of a 
family of functions, so that for 
constraint functions gj : 


KS(gj) 


1 

(1//>)LN(2 exp( P 

j-1 


gj)) 


(A 8 ) 


The KS function has the property of 


max (g-j) < KS(gj) < max(gj) + LN(l)/p 

(A9) 


Case 1: F easible Initial Design. 

Initial design data: (A17) 

X s - 10 in. for all b and h values 
in each finite element 
X c =0.01 
X s l = 3.0 in. 

X C L = 0.01 

W t = 5064.17 lb 

Optimal design data: 

*s 


X s u = 36.0 in. 

X c u = 0.06 

(A18) 

X C 


The effectiveness of that function in 
structural optimization has been reported 
many times, for example, in ref. 6 . In 
this application the KS-based cumulative 
constraints are also used in eq.A 6 for 
constraints g 6 and gj to represent 
dynamic displacements due to dynamic 
force F at several discrete time 
intervals. 

The design variable and behavior 
variable vectors for the structure and 
control subsystems are: 


3.1519 0.06 

14.177 

3.1133 

12.394 

3.0000 

10.141 

3.0000 

6.9576 

3.0000 

3.3719 

W t « 1509.17 lb* 


X s 


(b,h)x,(b,h)2,...,<b,h) 5 ) T 

(A10) 


*(4% violation of static stress 
constraint) 


w l 2 » w 2 2 » w 3 2 . *7- • • • / 

*36- 

Ws t 9l / • • • 9 95 ) 

(All) 

2 n 

c ) - 

(A12) 

u i 

n>c- w c , g 6 , g 7 ) T 

(A13) 


Case 2: Infeasible Initial Design. 

Initial design data: (A19) 

X s = 5 in. for all b and h values in 
each finite element 
X c =0.01 

(Lower and upper bounds are same as 
in Case 1) 


where the control subsystem weight is 

approximated as a function of the control Optimal design data: 
effort by an analytical expression: 

X s 


q 

W c = S K c I Ui |V2 ( A14 ) 


i=l 



3.1480 




14.255 

The numerical 

data for the example 

is: 

3.0123 




12.451 

E 

10.5 * 10 6 psi 

(A15) 

3.0430 

P s 

0.1 slug/ in 3 


10.115 

u 

0.3 


3.0000 

L 

500 in 


6.9571 

Tl 

1000 lb 


3.0000 

t 2 

5000 lb 


3.3197 

t 3 

1000 lb 






W t = 1502.48 lb 

The excitation 

force is given as a 

ramp 


function of time such that 



F (t) = k * 

t 

(A16) 


where k = 1000 

lb/sec. and the time 

1 t. 


varies from 0 

to 2 sec. The entire 




optimization problem is limited to the 
two second time period for the dynamic 
response. 


(A20) 

x c 

0.06 
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TABLE 1 Structures - Controls Coupling Data and Design Variables. 

Structures Subsystem 
(Equations: 1,2) 


Y s - output from Structures input to Controls: 

$ - matrix of natural vibration modes, 

w - array of the natural vibration frequencies. 


Y s - output from Structures to outside: 

W s - structures weight 

9l“95 “ values of constraints on static stress, static 
displacements and natural frequencies 


X s - design variables directly influencing structure: 

w,h - width and height of the beam rectangular cross-section 


Controls Subsystem 
(Equations: 3 through 13) 


Y c - output from Controls input to Structures: 

m c - control subsystem mass participating in the 
beam vibration 


Y c - output from Controls to outside: 

W c - control subsystem weight 

96 '97 “ values of dynamic displacement constraints 

X c - design variables directly influencing Controls: 

c - ratio of critical damping factor to frequency 


TABLE 2 Control Level: Low 

Finite GSE % 

Difference Difference 


DYs/DXc 

0 . 63960E+01 

0.63887E+01 

0.11 

0 . 20304E+03 

0 . 2018 4E+03 

0.59 

0 . 58682E-02 

0.85049E-02 

0.74 

0.26325E-01 

0. 26024E-01 

1.14 

0.X5871E-03 

0.15833E-03 

0.24 

0 . 32208E-03 

0 . 32118E-03 

0.28 

0 . 29355E-01 

0.29358E-01 

0.01 

0.31789E-01 

0. 31317E-01 

1.49 

0.24835E-03 

Q.24771E-03 

0.26 

-.31898E-03 

-.32132E-03 

0.73 

0.56326E-01 

0.56031E-01 

0.52 

42717E-01 

-.42761E-01 

0.10 

0 . 27978E-03 

0 . 277 13E-03 

0.95 

-. 11067E-02 

10984E-02 

0.76 

0.84043E-01 

0.83459E-01 

0.69 

- . 16504E+00 

”• 16436E+00 

0.41 

0.26853E-03 

0. 267 34E-03 

0.45 

- . 11968E-02 

1194 0E-02 

0.23 

0 . 11007E+00 

0. 10937E+00 

0.64 

- . 2668 3E+00 

-.26588E+00 

0.36 

0 . 254 56E-03 

0.25362E-03 

0.37 

— . 88010E-03 

87980E-03 

0.03 

0 . 00000E+00 

0. OOOOOE+OO 

0.00 

0. OOOOOE+OO 

0. 00000E+0O 

0.00 

0. OOOOOE+OO 

0. 00000E+00 

0.00 

- . 3 1988E-01 

31760E-01 

0.71 

- . 12755E+00 

12707E+00 

0.38 

0 . OOOOOE+OO 

0. OOOOOE+OO 

0.00 


PYg/PX<? 


— . 54 166E+00 

54 149E+00 

0.03 

-. 20907E+03 

20905E+03 

0.01 

0 . 88215E-01 

0.89733E-01 

1.69 

0. 13431E+00 

0. 13612E+00 

1.33 


TABLE 3 

Control Level: Moderate 

Finite 

GSE % 

Difference 

Difference 


PYs /PXC 


0 . 99818E+00 

0. 10989E+00 

9.17 

0. 34587E+02 

0 . 38611E+00 

10.42 

0 . 62088E-03 

0. 14330E-02 

56.67 

— . 56674E-01 

0. 50202E-02 

108.86 

0.82849E-04 

0 . 26614E-04 

67.88 

0 . 51425E-02 

0 . 62056E-04 

98.79 

0 . 46690E-02 

0.49235E-02 

5.17 

0 . 31988E-01 

0 . 60925E-02 

80.95 

- . 38805E-06 

0 . 41292E-04 

100.94 

- . 33497E-02 

60622E-04 

98.19 

0.83447E-02 

0 . 934 19E-02 

10.67 

- . 24339E-01 

— . 78102E-02 

67.91 

0 . 6907 3E-04 

0 . 45602E-04 

33.98 

0 . 21389E-02 

- . 20424E-03 

109.55 

0 . 12 318E-01 

0.13824E-01 

10.89 

— . ?3010E-01 

- . 29951E-01 

23.18 

0. 15522E-04 

0 . 43346E-04 

64.19 

- . 20629E-02 

- . 21153E-03 

89.75 

0. 15895E-01 

0 . 17997E-01 

11.68 

- . 19669E-01 

- . 47254E-01 

58.38 

0.51999E-04 

0 . 40773E-04 

21.59 

0. 12511E-02 

-.14467E-03 

111.56 

0. 00000E+00 

0. 00000E+00 

0.00 

0. OOOOOE+OO 

0. OOOOOE+OO 

0.00 

0. 00000E+00 

0. 00000E+00 

0.00 

- . 4967 IE-02 

- . 52758E-02 

5.80 

— . 21060E-01 

“ . 23602E-01 

10.77 

0. 00000E+00 

0. 00000E+00 

0.00 


P Y<? /P X <?_ 


-.77337E-01 

- . 81238E-01 

4.80 

- . 29844E+02 

- . 31440E+02 

5.08 

0.66012E-01 

0 . 87562E-01 

24.61 

0.79026E-01 

0 . 11993E+00 

34.11 
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TABLE 4 Control Level: High 


Finite 

GSE 

% 

Difference 


Difference 


ma/BXc 


0 . 78328E+01 

0. 78276E+01 

0.07 

0.29663E+03 

0 . 29595E+03 

0.23 

0.10070E-01 

0.10081E-01 

0.11 

0.38991E-01 

0. 38653E-01 

0.87 

0 . 18704E-03 

0. 18706E-03 

0.01 

0. 44703E-03 

0.47787E-03 

6.45 

0 . 34322E-01 

0. 34486E-01 

0.48 

0. 46889E-01 

0.47030E— 01 

0.30 

0 . 28638E-03 

0. 28819E-03 

0.63 

-. 43656E-03 

- . 46235E-03 

5.58 

0.64969E-01 

0. 65184E-01 

0.33 

58313E-01 

58194E-01 

0.20 

0 . 3 1665E-03 

0.31516E-03 

0.47 

15514E-02 

-. 15342E-02 

1.11 

0.96162E-01 

0 . 96032E-01 

0.13 

— . 22199E+00 

- • 22184E+00 

0.07 

0 . 29802E-03 

0.29639E-03 

0.55 

- . 15196E-02 

-.15293E-02 

0.64 

0.12438E+00 

0. 12448E+00 

0.08 

- . 34272E+00 

— . 34279E+00 

0.02 

0.28095E-03 

0 . 27663E-03 

1.54 

98255E-03 

— . 97852E-03 

0.41 

0 . OOOOOE+OO 

O.OOOOOE+OO 

0.00 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.00 

0. OOOOOE+OO 

O.OOOOOE+OO 

0.00 

36558E-01 

— . 36709E— 01 

0.41 

— . 17767E+00 

— . 17709E+00 

0.33 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.00 


■PYC/PXg. 


- . 52992E+00 

- . 52977E+00 

0.03 

-.20453E+03 

— . 20448E+03 

0.03 

— . 84 192E-01 

- . 7897QE-Q1 

6.20 

- . 16570E+00 

— ♦ 16844E+00 

1.63 



Figure 1.- Cantilever beam with actuators. 



Figure 3.- Histogram of optimization started 
from a feasible design. 


Weight 

Ob) 



Iteration number 


Figure 4.- Histogram of optimization started 
from infeasible design. 



Figure 2.- Coupled, structure-control system. 
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